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We discuss numerical results of relativistic magnctohydrodynamic (MHD) jet formation models. We first reviews some 
examples of stationary state solutions treating the coUimation and acceleration process of relativistic MHD jets. We 
provide an a posteriori check for the MHD condition in highly magnetized flows, namely the comparison of particle 
density to Goldreich-Julian density. Using the jet dynamical parameters calculated from the MHD model we show the 
rest-frame thermal X-ray spectra of the jet, from which we derive the overall spectrum taking into account a variation 
of Doppler boosting and Doppler shift of emission lines along the outflow. Finally, we present preliminary results of 
relativistic MHD simulations of jet formation demonstrating the acceleration of a low velocity (0.01 c) disk wind to a 
coUimated high velocity (0.8 c). 



1. Introduction 

Astrophysical jets are commonly believed to consist 
mainly of coUimated disk winds, launched, acceler- 
ated, collimated by magnetic forces (see [1-5, 7-9]). 
The jet is initiated as a slow wind by a process which 
is not yet completely understood, in particular its 
time-dependent character. Most probably, some disk 
instability is responsible for adjusting the mass flux 
from disk into an outflow in the direction perpen- 
dicular to the disk surface. The disk wind is first 
launched magneto-centrifugally and further acceler- 
ated and (self-) collimated into a narrow beam by 
Lorentz forces. It should be noted, however, that cer- 
tain important properties of relativistic jets are still 
purely known. We don't know whether relativistic 
jets are launched as leptonic or hadronic beams. We 
do not know much about the intrinsic jet velocities. 
We also don't know the magnetic field structure of 
the jet or of its source. 

Traditionally, the theoretical task to settle these 
issues can be divided into five fundamental questions: 

• The question of how to coUimate and accelerate a 
disk/stellar wind into a jet (jet formation)? 

• The problem of how to eject the disk/stellar ma- 
terial into the wind (jet launching)? 

• The question of the accretion disk structure and 
evolution? 

• The question of the origin of the governing mag- 
netic field? 

• The question of the jet propagation, its interac- 



tion with ambient medium, its stability and the 
radiative signatures involved. 

In the end, it would be desirable to obtain fully self- 
consistent solutions of relativistic MHD jets - how- 
ever, such solutions have not yet been obtained in the 
stationary case. However, the progress in numerical 
methods and computer power now allows for time- 
dependent simulations of relativistic jet formation. 

In this work we concentrate on the first prob- 
lem, the structure and dynamics of the jet forma- 
tion region close to the source. We present stationary 
state and time-dependent numerical results of rela- 
tivistic jets as collimating disk winds. For alternative 
or complementary viewpoints (e.g. electromagnetic 
jet formation, Blandford-Znajek mechanism) we re- 
fer to the literature. 

2. Concepts of relativistic MHD jets 

Considerations of stationary MHD flows have re- 
vealed that relativistic jets must be strongly mag- 
netized [3, 10-12]. In that case, the available mag- 
netic energy can be transfered into a small amount of 
mass with high kinetic energy. On the other hand, a 
very strong magnetization may be in conflict with the 
MHD assumption lacking a sufficient large amount of 
electric charges which are needed to drive the electric 
current system. 

Theoretical modeling of jet formation requires to 
solve the governing MHD equations. However, due to 
the complexity of MHD and the astrophysical bound- 
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ary conditions indicated, a completely self-consistent 

MHD solution for the jet formation process being 
compatible with all the generic features (MHD self- 
collimation, accretion-ejection mechanism, magnetic 
field generation, spatial and time scales etc.) does 
not yet exist. 

We believe that the jet formation region may 
fairly well be approximated as axisymmetric. Non- 
axisymmetric distortions may actually hinder the 
formation of powerful jets as probably demonstrated 
by the existence of many strongly magnetized, 
rapidly rotating accretion disk systems which do not 
have jets (e.g. cataclysmic variables or most pulsars). 

Compared to the non-relativistic case, in rela- 
tivistic MHD some specific features arise affecting 
both the model approach and the physical mech- 
anism of jet formation. We first mention the light 
cylinder of the magnetosphere. If ftp is interpreted 
as angular velocity of magnetic field lines, at the light 
cylinder radius _Rl = c/D,-p of that field line its ro- 
tational velocity equals the speed of light. In general 
relativity, frame dragging adds another, inner, light 
"cylinder" to the system of MHD equations. For van- 
ishing mass density, the location of the MHD Alfven 
surface approaches the light surface. Another rela- 
tivistic MHD feature is that the electric field can- 
not be neglected. The poloidal electric field com- 
ponent E± (which is perpendicular to the poloidal 
magnetic field) reaches considerable strength, E± = 
{R/RL,)Bp. The electric field pressure gradient and 
tension adds to the jet coUimation or de-collimation. 

We like to note that neither the light surface 
nor the relativistic Alfven surface obeys a self-similar 
structure (see e.g. [11]). Thus, steady-state self- 
similar relativistic MHD models while having the 
advantage of treating the matter-field interrelation 
self-consistcntly (compared e.g. to force-free models 
discusscid below) may not necessarily take into ac- 
count all relativistic properties of the flow. 

3. Stationary relativistic MHD jets 

In stationary limit global solutions for the relativis- 
tic jet structure can be derived (see e.g. [11, 13-15]) 
on spatial scales and with a resolution which can- 
not (yet) be reached by time-dependent simulations 
[16, 17, 30]. Time-dependent MHD codes usually 
have also difficulties with low plasma-/3, and/or high 
magnetization - exactly the conditions which gener- 
ate highly relativistic jets. On the other hand, jets 
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Fig. 1. Solutions of the stationary state axisymmetric MHD 
equations in Kerr metric. Axisymmetric poloidal magnetic 
field structure (top, sec [14]), cylindrical coordinates normal- 
ized to gravitational radii. The inner boundary coincides with 
the inner light surface. Maximum poloidal velocity along a col- 
limating flux surface Mp,oo as a function of flow magnetization 
a {bottom, see [22]), both in Kerr metric. 

are non-steady phenomena as clearly seen by their 
knotty pattern. 

Example solutions of stationary state (general) 
relativistic MHD are shown in Fig. 1. The figure 
shows the global axisymmetric field structure of the 
jet collimation region considering the (force-free) lo- 
cal cross- field force-balance (see [14]). This is a truly 
2.5-D solution, taking into account special relativis- 
tic and general relativistic features (light surfaces, 
gravitational redshift, frame dragging). The asymp- 
totic structure (for large z) matches perfectly with 
the analytical 1-D solution. The electric current is 
conserved along the field lines. The global field solu- 
tion (inside and outside the light cylinder) is deter- 
mined by matching by the regularity condition (along 
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Fig. 2. Solutions of the stationary state relativistic MHD 
wind equation for a coUimation field structure. Shown is the 
ratio of particle density n to Goldreich-Julian space charge 
density uqj versus cylindrical radius along the flux surface 
(see [19]) for examples of highly magnetized flows (top) with 
a = 1000 and cr = 5000 (bottom). The physical branch of the 
wind solution is the one which starts with high value at small 
radius and decreases with radius. 

the light cyhnder). EssentiaUy, the field structure is 
coUimating rapidly. Similar solutions were obtained 
for Minkowski case, however assuming differentially 
rotating foot points of field hnes [18]. In that case, 
the light surface has to be determined iteratively as 
its location and shape is not known a priori since 
both depend on the magnetic flux distribution. 

Figure 1 also shows an example for the flow dy- 
namics along collimating field lines, considering the 
force-balance along the flow and taking into account 
magnetic and inertial forces and also gas pressure. 
Shown is the maximum poloidal velocity along a 
field line (taken at a certain radius) as a function 
of the flow magnetization (see [19]). Essentially, the 
velocity-magnetization relation follows a variation 
of the well-known Michel scaling [10]. This exam- 
ple solution was calculated for hypothetical jets in 
GRBs. For high magnetization, bulk Lorentz factors 
of 1000 were obtained. The transition between al- 
most conical field lines (split-monopole) close to the 
source and asymptotic cylindrical field lines leads to 
an interesting behavior. As these flux surfaces ac- 
tually re-coUimate (see the separation between the 
flux surfaces), for high magnetization we do not flnd 
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Fig. 3. Thermal X-ray spectra from relativistic MHD jets. 
Shown are signatures of the underlying magnetohydrodynamic 
structure - the poloidal jet velocity Up = yUp along the flux 
surface (top) and the gas temperature T (bottom), see [22, 24]. 
Distance x is normalized to the gravitational radius. 

stationary solutions for each radius beyond the re- 
collimation distance. This indicates that the flow is 
highly time- variable, giving rise to internally differ- 
ential motions and thus internal shocks which may 
trigger the prompt emission of the GRB (see [19]). 

4. Relativistic jets &: MHD assumption 

An essential issue for relativistic MHD models is the 
question whether the MHD limit actually holds for 
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Fig. 4. Thermal X-ray spectra from relativistic MHD jets. 
Shown are example spectra derived for the accelerating hot 
gas. The top panel shows a hot (lO'^ K), slow gas component, 
the lower panel a cooler (10^ K), but faster gas component 
[24]. The different lines show the spectra for diflferent viewing 
angles (-40, -20, 0, +20, +40 degrees). 

the parameter space investigated. A small mass load 
implies a vanishing amount of charge carries. Elec- 
tric currents may not longer be established. The va- 
lidity of the MHD condition can be estimated by 
comparing the particle density n in respect to the 
Goldreich-Julian density ncj (see [19, 20]; see [21] for 
another approach). Figure 2 demonstrates our a pos- 
teriori check for the MHD condition for highly mag- 



netized flows. For large radii the fraction (n/ncj) 

becomes constant value ^ 100 >> 1, telling us that 
our MHD solution for this choice of parameters is 
consistent with the MHD assumption. The solutions 
shown in Fig. 2 were derived for coUimating field lines 
with slightly diverging flux, BpV^ = ^{r) ^ r'^ with 
g ~ 0.1 > 0. For a more diverging field structure 
(lower q), a radius will exist beyond which the MHD 
condition will break down, {n/ncj) — 1- 

5. Thermal X-rays from relativistic jets 

Accretion disks around stellar mass black holes may 
reach temperatures up to 10^ K close to their in- 
ner edge. Assuming that this disk material becomes 
ejected into a disk wind and accelerated to a jet, 
each volume element of matter may radiate thermal 
Bremsstrahlung in the 1-10 keV X-ray band. In case 
of optically thin emission, emission lines will be vis- 
ible in the corresponding spectrum. Relativistic mo- 
tion of each volume element will boost and Dopplcr- 
shift its rest-frame spectrum. As the jet coUimates, 
the line-of-sight angle of the moving elements is 
changing, leading to a variation in the Doppler factor 
along the jet. 

Figures 3 and 4 show example steps of our 
derivation (see [24] for details). In the first step we 
calculate the MHD dynamics of a coUimating jet by 
solving the " hot" wind equation in Kerr metric from 
about 5 to 500 gravitational radii [22]. This provides 
temperature, density, velocity etc. along the flow. In 
the second step we calculate the rest-frame spectrum 
of each of the fluid elements along the jet, each hav- 
ing different temperature, density (Fig. 4). For ther- 
mal Bremsstrahlung the continuum is dominant for 
temperatures > 10^ K. A number of lines appear for 
lower temperatures/densities (e.g. lines of O, N (0.5- 
0.9 KeV), Ne, Fe, Mg (1-4 KeV) or FcXXV, FeXXVI 
emission (6.6-7.0 KeV). In the last step we combine 
the local rest-frame spectra of all fluid elements, each 
moving with a different velocity in respect to the line 
of sight and, thus, affected by a different Doppler fac- 
tor. For our example of a micro-quasar we find a total 
X-ray luminosity in the rest frame of 4 x 10^^ ergs/s, 
respectively 6 x 10"^^ ergs/s beamed and boosted and 
seen from along the jet axis, or 2 x 10^'^ ergs/s as seen 
from 20 degree inclination. This is much less than 
the kinematic jet luminosity of Lkin — 10^^ ergs/s 
(assuming mass loss rates of 10~^° M0/yr). 



November 20, 2008 22:1 



WSPC - Proceedings Trim Size: llin x 8.5in fendt 



5 



6. Simulations of relativistic MHD jets 

Finally we present preliminary results of numerical 
simulations considering the formation of relativistic 
MHD jets. The overall goal is to test whether the 
paradigm of MHD self-collimation known for New- 
tonian jets [5, 25, 26] also holds in the relativistic 
case. 

A Keplerian disk is taken as a fixed-in-time 
boundary condition with certain magnetic flux pro- 
file and mass loss rate from disk into corona. This ap- 
proach allows for long-term simulations of the evolv- 
ing jet which is the main advantage of not considering 
the disk evolution in the simulation. For simulations 
taking into account the disk structure we refer to the 
literature [27-30]. 

Our ideal MHD simulations start from the ini- 
tial state of a hydrostatic corona within a potential 
magnetic field (see [31]). We assume a polytropic gas 
with polytropic index 7 = 4/3. Turbulent magnetic 
difl[usivity [26] or the variation of the magnetic flux 
and mass flow profiles across the jet [32] is not yet 
taken into account. Mass is "injected" from the disk 
surface into the jet with low velocity of about 0.01 
times the local Keplerian speed VK{r)- 

We use the new MHD code PLUTO [33] in the 
special relativity option adding a Newtonian descrip- 
tion of gravity. The grid is non-equidistant with high 
resolution close to the jet axis and the disk surface, 
spanning 40 x 40 inner disk radii Rin- The leading 
dynamical parameters for the simulations presented 
here are a plasma-/? = 0.1, VK{r = Rin) = 0.4c. 

Figure 5 shows our preliminary results. Shown 
is the axial component of the velocity Vz {r, z) jc for 
several time steps. The low velocity material injected 
from the disk surface becomes accelerated to about 
0.8 c. We observe two jet components - a coUimated 
high speed beam emerging from the inner disk and an 
outer low velocity disk wind. However, while the sim- 
ulation reaches a quasi stationary state close to the 
jet axis, the outer regions still dynamically evolve. 
This is natural as 100 rotations at the inner disk ra- 
dius correspond to 0.5 rotations at the outer disk. 

The whole outflow is still sub-Alfvenic with low 
toroidal magnetic field strength. The degree of col- 
limation as measured by the directed mass fiux is 
about 5-10 degrees. Since the flow is sub-Alfvcnic we 
assume coUimation by the surrounding low velocity 
wind. 
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